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Abstract In this study, we consider a flat Friedmann- 
Robertson- Walker (FRW) universe in the context of 
Palatini f(R) theory of gravity. Using the dynamical 
equivalence between f(R) gravity and scalar-tensor the- 
ories, we construct a point Lagrangian in the flat FRW 
spacetime. Applying Noether gauge symmetry approach 
for this f(R) Lagrangian we find out the form of f(R) 
and the exact solution for cosmic scale factor. It is 
shown that the resulting form of f(R) yield a power- 
law expansion for the scale factor of the universe. 

Keywords f(R) gravity; Palatini formalism; Noether 
gauge symmetry 
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1 Introduction 

It is well-known that the Einstein-Hilbert action pro- 
duces general relativity (GR) which leads mostly to 
observational success. Furthermore the modified grav- 
ity theories issued by more general gravitational ac- 
tions could explain the observation al facts. Obser- 
vations of the superno vae Type la ( Riess et al. 19981 : 
Perlmutter et al. 19991) and the cosm ic microwave back-| 
ground radiation (jNetterfield 20021 ) indicate that cur- 
rent expansion of the universe is accelerating, on the 
contrary to Friedmann-Robertson- Walker (FRW) so- 
lution of GR. To explain such an accelerated expan- 
sion, in the framework of GR, many authors intro- 
duced mysterious cosmic fluid, the so-called dark energy 
which can be described by the cosmological cons tant 
( Copeland et al. 20061 iDurrer and Maartens 20081) . 
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On the other hand, to get a solution of this prob- 
lem, some modifications in GR theory have been pro- 
posed. The f(R) gravity th eory is a particular class of 
modified gravity theories (ISotiriou and Faraoni 20101: 
Noiiri and Odintsov 20o"7tlDe Felice and Tsuiikawa 2010[ ).| 
The simplest form of this theory can be constructed 
by replacing the Ricci scalar R with an arbitrary func- 
tion f(R) in the Einstein-Hilbert Lagrangian (EHL). In 
recent literature some different forms of f(R) gravity 
have been proposed, and discussed in different con- 
texts. For example, it has been shown that early- 
time inflation and current cosmic acceleration may take 
place by adding ne gative and positive powers of curva- 
ture int o the EHL dNoiiri and Odintsov 20031 ). Carroll 
et al. (jCarroll et al. 2004 I ) have also proposed that 
a small correction to the EHL by adding an inverse 
term of R would lead to cosmic speed-up which orig- 
inates from purely gravitational effects. Similar mod- 
ification s of GR have also been proposed to drive in- 
flation ( Starobinskv 198oh . On the other hand, it is 
worth noting that the main deficiency of such theo- 
ries is tha t they are se r iously constrained by solar sys- 
tem tests jOlmo 2005at IChiba et al. 20071) . A number 
of viable f(R) models that can satisfy both cosmologi- 
cal and loc al gravity constraints h ave been proposed 
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in Refs. (lAmendola et al. 2007 
Cognola et al. 20081 ). 

There are two kinds of Noether symmetry approach 
for cosmological studies in the literature : The first 
one is the so-called Noether symmetry approach in 
which the Lie derivative of a giyen Lagrangian vanishes, 



( Capozziello and de Ritis 19931; de Ritis et al. 1990 ; Demianski et al. 1 



Sanval and Modak 20011 : ICamci and Kucukakca 2007 ; 
Souza and Kremer 20081 ) and the second one is the 
so-called Noether gauge symmetry (N GS) approach 
( Jamil et al. 20 111: iHussain et al. 201 lh . The latter is 
a generalization of the former Noether symmetry ap- 
proach in the sense that Noether symmetry equation 
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includes a gauge term. Taking into account a gauge 
term in Noether symmetry equation gives a more gen- 
eral definition of the Noether symmetry: this is the 
so-called Noether gauge symmetry (NGS) approach. 
Thus one may expect extra (more than one) symmetry 
generators from this definition. 

We note that the Noether symmetry approach 
without gauge term allows one to choose the poten- 
tial dynamically in the scalar-tensor gravity theory 
( Sanval et al. 20031 ) . and very recently, the explicit 
form of the f unction f(R) (|Capozziello and de Felice 2005 



Vakili 20081) . For this form of f(R), cosmological so- 



lutions of FRW metric can describe the accelerated 
period of the universe. Applying this Noether sym- 
metry approach, the spherically symmetric solutions 
in f(R) theories of gr avity have been also found in 
ICapozzillo et al. 20071). Using the definition of NGS 
some authors ( Jamil et al. 20 lit Hussain et al. 2011 ) 
have discussed f{R) and /(i?)-Tachyon model in the 
metric formalism, where they conjectured that the ap- 
plication of NGS to generic f{R) Lagrangian results 
in zero gauge function. The Palatini formalism has 
not been considered by using the NGS approach in the 
literature yet. In a rece nt work Roshan and Shojai 
(jRoshan and Shoiai 20081 ) studied Palatini f(R) cos- 
mology in flat FRW spacetimes following Noether sym- 
metry approach without gauge term for the matter 
dominated universe. They found out the form of f(R) 
as power-law and exact solutions for cosmic scale factor. 
In this study we have generalized their result consid- 
ering NGS approach within the scope of the Palatini 
f(R) gravity. 

This paper is organized as follows. In the following 
section, the Palatini formalism is briefly reviewed. In 
section [3j we discuss the Noether symmetry approach 
with/without a gauge term for the Palatini f(R) grav- 
ity. In section 21 we search the cosmological solutions 
by using the obtained forms of f(R)- Finally, in Section 
[5J we conclude with a brief summary. 



2 The f{1Z) Gravity and The Field Equations 

In four dimensions, the action of the Palatini f(lZ) grav- 
ity theory with matter is written as 



A 



d i x^[f{'R.)+2KL m (g ab ,^)}. 



(1) 



Here / (TV) is a differentiable function of the Ricci scalar 
TZ = g ab Tt a b(T), lZ ab (T) is the Ricci tensor of any 
independent torsionless connection T independent of 
g a b and L m is the matter Lagrangian and ip repre- 
sents collectively the matter fields. The matter La- 
grangian is chosen as L m = —p m oa~ 3 for matter domi- 
nated cosmology. In general there are two approach in 



order to derive the dynamical field equations of mo- 
tion in f(R) gravity. The first on e is Palatini ap- 
proach Flanagan 20041 loimo 2005bl iFav et al. 2007 



Baojiu et al. 20071 ) in which the metric and connection 



are considered as independent quantities, and the ac- 
tion is varied with respect to both of them. The sec- 
ond approach is the metric formalism in which action 
is varied with respect to metric tensor. The field equa- 
tions in the metric formalism are fourth-order differ- 
ential equations, while for Palatini formalism they are 
'Isecond-order. If f(R) is linear in R, the two approaches 
lead to the same equation. In this study we will use the 
Palatini formalism. Variation of the action (JlJ with re- 
spect to metric t ensor yields following field equations 
( Brax et al. 20081) 



faTZ-ab - -^fgab ~ K-T abl 



(2) 



and the variation of the action ([T]) with respect to the 
connection gives 



a b \ 



o, 



(3) 



where fa — df/dlZ and T ab is the usual stress-energy 
tensor of the matter. Also, the trace of Eq. is 



fan - V = kt. 



(4) 



The f(TV) gravity theory in Palatini formalism is equiv- 
alent to lobd = —3/2 a Brans-Dicke (BD) theory. In 
order to construct a canonically effective point-like La- 
grangian, we have to use the dynamical equivalence be- 
tween P alatini f(R) f o rmalism and the BD theory of 
gravity |Sotiriou 20061 ICapozzillo et al. 201ll ). There- 
fore the action ((TJ) can be written as follows 



A 



d i x^g-^R+—g ab d^ 



%<f>-U(4>)+2 K L m ), (5) 



where = fa, U(<f>) = 0x0) ~ f{x{<t>)), R = x{<t>) and 
R is a curvature scalar constructed from the Levi-Civita 
connection of metric tensor. We note that this is the 
action of the BD theory with the BD parameter lobd = 
—3/2. It was also shown that in the metric formalism 
the f(R) gravity is equivalent to th e BD theory with 
ojg £) = (jSotiriou and Faraoni 20101 ) . 

In this study, we consider a matter-dominated 
model in a spatially flat FRW universe with signature 
(+,+,+,—). Using the action ([5]) with the scalar curva- 
ture for FRW metric in the form R = 6 (a/a + a 2 /a 2 ) 



a point-like L agrangian 



and taking L m 

takes such a form ( Roshan and Shoiai 20081 ) 



3a • 

L = 6a 2 a<j> + 60ad 2 + — <j) 2 + a 3 C/(0) + 2np m0 , (6) 
20 
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where a(t) is the cosmic scale factor, p m Q is an integra- 
tion constant associated with matter content and the 
dot indicates the derivation with respect to cosmic time 
t. It is found that the energy function El associated 
with the Lagrangian ((5]) vanishes, i.e. 



E, 



a atfi l<f> U Kpmo 



acb 4(j> 2 6</> 3a J 



0, 



(7) 



which is known as the modified Friedmann equation. 
The equations of motion can be obtained by varying the 
Lagrangian © with respect to a and <b, respectively, as 
follows 



2a a 2 (f> 2a<j) 30 2 



U 



(j) a<\> ±<t> 2 20 



0. 



<\> 3a<t> <j) 2 



a 2 2(f> 2a<p 



(8) 



(9) 



For Noether symmetry without gauge term (i.e. B — 0) 
it is required that r = 0, and thus the above equation 
reduces to the form £ X [i]L = which is the existence 
condition for Noether symmetry without gauge term. 
The si gnificance of N GS is clear from the following the- 
orem jNoether 19181 ). 

Theorem: If X is the Noether symmetry correspond- 
ing to the Lagrangian L(t, a, <j), a, </>), then 

F)T Fit 

I = tL + (a -to)Z- + 09-r^-fl (15) 
oa 

is a first integral or a conserved quantity associated 
with X. Now we seek the condition in order that the 
Lagrangian ((6]) would admit NGS. 

For the flat FRW metric, the NGS condition (0 
yields the following set of equations 



0, 



(16) 



3 The Noether symmetry approach 

For most of studies in the context of both f(R) grav- 
ity and the scalar-tensor theory, it has been used def- 



( CaDOzziello anc 


de Ritis 1993: Kamilya et al. 2004: 


Vakili 2008h . ( 


Roshan and Shoiai 2008). Recentlv. 



some work on NGS approach for FRW metric ha s 



been appeared (jJamil et al. 2011tlHussain et al. 20111 ). 



Noether gauge symmet ry is defined as follows. Let us 
consider a vector field (jlbragimov 1 999) 



X — t ^ + a® +/3^ 
dt da dch 



(10) 



where r, a and (3 are depend on t, a, and <j>. Here, t is the 
independent variable, a(t) and <f>(t) are the dependent 
variables. The first prolongation of the above vector 
field is given by 



da d<f> 

in which 

at = D t a - aD t T, j3 t = D t f3 - (j>D t T, 



(11) 



(12) 



where D t is the operator of total differentiation with 
respect to t 



d d ■ d 
Dt= d-t +a d-a + ^W 



(13) 



The vector field X is a NGS of a Lagrangian L(t, a, </>, d, < 
if there exists a gauge function, B(t, a, 4>), such that 



X [1] L + LD t (r) =D t B. 



(14) 



Qa{2<t>a t + af3 t ) -B a = 0, (17) 
3a 2 

— {2</>at + ap t ) - Bj, = 0, (18) 



2a + a (a„ + - n) + 2^ + — f3 a = 0, (19) 



cf ) a + al3 + a ( f>{2a a -T t ) + a 2 l3 a = Q, (20) 



3a- -/3 + a(2/3 -r t ) + 40a = 0, (21) 



a 2 (3a + ar t )U + a 3 U'/3 + 2np m0 T t -B t = 0. (22) 

After some algebraic calculations, the solutions of the 
above set of differential equations (| 17 1) - ([22]) for a,f3, 
t, B and the potential U (</>) are obtained as follows: 



, 2 ,n-4±L 

a = c\a + cia\a <p) 2e + 1 , 

p = -(2e-i)<i> 



3C1 . 2 ,s_J±i. 



(23) 
(24) 



3ci(£ + l) 6 Cl « Pm0 (l+l) , 
— - — t + c 3 , B = — 1 + c 4 (25) 



C/(0) = A^ — , 



(26) 
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where Cj, £ and A are constants and I ^ —1/2, 2. Since 
B = and r = for Noether symmetry without gauge 
term, the parameters ci,C3 and C4 should vanish. Thi s 
case has been studied in Ref. (|Roshan and Shoiai 20081 ). 
Taking (f> = tp 2 in the Lagrangian © and doing the 
above calculations again, the obt ained form of a and 



it is possible to solve f(TV), where 1Z = xi'P) an( i 
fn- Putting the potential (|26[) into Eq. (|35[) yields 



fcK-f = \{M)-*=i. 



(36) 



Thus, it is straightforward to get the following two so- 



generalizes the ones found in Ref.( |Roshan and Shoiai 2008h | ltions from E ^ ® for fCQ 
(which are represented as A and B in this reference). 

It is seen from (|23)) - (l25l that the Lagrangian (|6|) ad- 
mits three Noether symmetry generators 



Xi = 



at' 



(27) W) = R Tl - XR 
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d 



X 2 = (a^)-(^) ( a— + {21- 1)0— 



(28) 



3(l+l)t 9 9 3(21 -!)</> d 
X3 = T^~9t + a d~a ~ 1-2 W (29) 

The corresponding Lie algebra has the following non- 
vanishing commutators: 



[Xi,X 3 ] 
[X 2 , X 3 ] 



3{£ + l) 



Xi, 



1-2 



Xo 



(t-2)(2t+iy 

The first integrals associated with Xi are 
3n 3 . 

ii = 6a"aq) + 6acf>a H -<fi — a U — 2Kp m o, 

2(f> 



h = 3(2^+ l)a(a 2 0)~(^> 4(« 2( i 

at 



(30) 
(31) 

(32) 
(33) 



(37) 



(38) 



A3 



-3/(21-1) 



l-2£ 



2JI+1) 



Ro is a con- 



where p{£) = 2{l + 1) 

stant and £ ^ -1,1/2. 

In the following subsections, considering f(lZ) given 
by ([37]) and (|38|) . we will search the exact solution for 
cosmic scale factor. For a and (3 given by (|23 ]) - (|24| and 
Eq. (|37|). it is appeared the singularities for the values 
of £ = —1/2 and I = — 1. For these values of £, after 
solving Eq. ([36| one has the same form with (|38| for 

W)- 

4.1 Case (i). f{TZ) = p{£)W^ 

Using the trace relation Q and <f> = fn we obtain 



</> = r(£)a 



2£-l 



(39) 



where r{£) = ^[^m)]^, £ ^ -1/4. Hence, 
the Noether first integrals (|33f and (|34|) can be written 

as 



3(2£ + lJ^^WTfl^ 



(40) 



h = 



3(£- 



£-2 



1) 3(4£+l)a d(a 2 



£-2 eft 
6(1 + l)Kp m0 
£-2 



t. 



(34) 



Here the constant parameter C4 is assumed to be zero 
in the gauge function B. We note that the first integral 
(|32|) is related with the energy function (|7|) , so that the 
first integral I\ vanishes. 



4 The Cosmological Solution 

As an inverse problem of finding f{lZ) Lagrangian, it is 
only required to give U (ip). Using the algebraic relation 



u(4) = M4>) - fbett)), 



(35) 



3(41+1) (21 +l)r(l) (2e+1) . 



£-2 



6np m0 (£+ 1) 
£-2 



t. 



(41) 



Now the Eq. (|4"0"|) can be used to find out the time de- 
pendence of the cosmic scale factor as 



a(t) = 



I 2 (l+l)r(l)-^+r 
3(2£+l) 2 



t + t 



(42) 



where to is a constant of integration. Considering a(t) 
in Eq. (|41[) , the following constraint equations arc found 



(43) 



Il(4£ + l)r(J0 1/(2m) - 18K PmQ (2£ + l) 3 = 0, 



I 2 t {U + l)r(£)( £+1 )/( 2£+1 ' + {£- 2){2£ + 1)J 3 = 0. (44) 



5 



It is explicitly seen here that if to = 0, then the first 
integral I3 is equal to zero. Also, using Eqs.(|7J), © and 
(151). wc have additional constraint relations as 



lj{2£-l)r{£) x l^+^ 



18(2^+1) 



Xr{£) 



-3/(21-1) _ 



0, 



(45) 



I$r(£) 



1/(2^+1) 



6(2^+1) 



\r(£)-W e -V +2K Pm0 \ =0. (46) 



Thus, the constraint Eqs. ([43)) and (|45|) give the fol- 
lowing simple relation 



Kp , n0 . _^±i) r(f) -3/ P <- 1 ». 



(47) 



The deceleration parameter, which is an important 
quantity in cosmology, is defined by q = —ad/ a 2 , where 
the positive sign of q indicates the standard decelerat- 
ing models whereas the negative sign corresponds to 
accelerating models and q = corresponds to expan- 
sion with constant velocity. It takes the following form 
in this model 



q = £. 



(48) 



The effective equation of stat e parameter defined by 
Wr.tf = -1 - = dCapozziello et al. 2003 : 



Allemandi et al. 2004a ; 
be obtained as 



Allemandi et al. 2004bf ) can 



2^-1 



,J eff 



(49) 



where H is Hubble parameter. Astrophysical data indi- 
cate that w lies in a very narrow strip close to w = — 1 . 
The case w = — 1 corresponds to the cosmological con- 
stant. For w < — 1 the phantom phase is observed, 
and for —1 < w < —1/3 the phase is described by 
quintessence. Thus, in the interval — 1 < £ < we have 
quintessence phase. If— 00 < £ < — 1, then the phan- 
tom phase occurs, where the universe is both expanding 
and accelerating. 



4.2 Case (ii). f(K) = R K - XR 



3/(2^-1) 



If i?o = 1, i-c = 1 from the relation <fi = fn = Rq, 
then the action is reduced to the Einstein-Hilbert action 
with cosmological constant ( f(lZ) = 1Z — A). Wc note 
here that Palatini and metric formalism are coincide. 
In this case, the Noether first integrals (|33|) and (f34|) 
can be written as 



h = 6(2£ + l)a~^a, 



6(4£ + T 



a a + 



6(1 + l)np m o 
1-2 



t. 



(50) 



(51) 



The modified Fricdmann equation ([TJ for this case re- 
duces to the form 



K-PmO 

3a 3 



0. 



From Eq. ([50]) the scale factor is solved as 
"I 2 (4^+l) 



o(f) 



21 + 1 
47+T 



6(2^+l) 2 



t + U 



(52) 



(53) 



Inserting the scale factor ([53)1 into Eq. ([51]) and Eq. 
((52)) one gets A = 0, £ = 1/2 and as a constraint equa- 
tions npmo = if/48 and t\ = /3/J2. Therefore, for this 
case the scale factor has the form 



a(t) 



h h 
—t + — 

8 h 



2/3 



(54) 



which is obtained from (|53[) when £ = 1/2 . Taking 
t = -1/4 in Eq.([5T)]). we have a(t) = exp(^) which 
gives p m o — and A = 2lf/3. Using these results for 
£ = — 1/4 in the Eq. ([51]) one can find I 3 = 0. 



5 Concluding remarks 

The Palatini approach consider the metric g a i, and 
the connection T^ c as independent field variables, but 
the spacetime metric g a b is only independent variable 
in metric formalism. The Palatini formalism can be 
seen as containing two independent metrics g a b and 
h-ah = fizgab rather that a metric and independent con- 
nection. In Palatini f(lZ) gravity the second metric h a b 
determine the geodesic structure with the connection 
r^ c which is the Levi-Civita connection of new metric 
h a b- In BD theory of gravity the second metric h a b is 
related to the non-minimal coupling of the BD scalar, 
i.e. (f) — f-R- In Palatini approach, f(1Z) gravity given 
by the action ([T]) is equivalent to a special BD theory 
with a scalar field potential. The action ([5]) is clearly 
that of a BD theory with BD parameter Wbd = — 3/2 
and a potential U (</>), which is considered in this pa- 
per. In the metric formalism, f(R) gravity is equiva- 
lent to the BD theory with wbq = (see the review 
paper ( Capozziello and Laurentis 20 111 )). An wbd = 
BD theory was originally studied for the aim of getting 
Yukawa correcti on to the Newt onian potential in the 
weak- field limit |0'Hanlon 1972h . 

In this study, we examined the matter dominated 
flat FRW universe by considering the Palatini f(lZ) 
formalism and by following the NGS approach which 
leads to explicit form for f(TZ). This approach is based 
on the search for Noether gauge symmetries which al- 
low one to find the form of f{TVj. We have obtained 
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two type Noether symmetric f(TZ), (|37| which yields 
case (i) and (|38|) which gives case(ii), where f(lZ) func- 
tions give rise to a power-law Lagrangian and EHL 
with cosmological constant, respectively. In case (i), 
it is found the power-law form of cosmic scale factor 
by (|4"2")l . For case (ii), the cosmic scale factor is ob- 
tained by ([54")) for £ = 1/2 and de Sitter solution for 
I = —1/4. We have presented the effective equation 
of state parameter for Palatini f(lZ) cosmology. In the 
first model, case(i), the expansion of Universe is accel- 
erating at the intervals — 1 < £ < 0, quintessence phase, 
and — oo < £ < —1, phantom phase. Thus, this model 
can provide a natural gravitational alternative for the 
dark energy without the necessity to introduce an ex- 
otic fluid with a negative equation of state parameter. 
We note here t hat while the gauge function turns out t o 



be zero in Refs. (| Jamil et al. 2011uHussain et al. 20111 ). 



but our analysis shows that it depends on the cosmic 
time. 
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